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Abstract 

We consider the problem of sound propagation in a wind. We note that the rays, as 
in the absence of a wind, are given by Fermat's principle and show how to map them to 
the trajectories of a charged particle moving in a magnetic field on a curved space. For the 
specific case of sound propagating in a stratified atmosphere with a small wind speed we 
show that the corresponding particle moves in a constant magnetic field on the hyperbolic 
plane. In this way we give a simple 'straightedge and compass' method to estimate the 
intensity of sound upwind and downwind. We construct Mach envelopes for moving sources. 
Finally, we relate the problem to that of finding null geodesies in a squashed anti-de Sitter 
spacetime and discuss the SO(3, 1) xR symmetry of the problem from this point of view. 

Keywords: 



1. Introduction 

In the New Scientist^ of 15 th April, 2009, a correspondent posed the following question 
regarding traffic noise 

I live a kilometre north of a busy motorway. When the wind is coming from 
the south the noise of the motorway is noticeably greater than when the wind is 
coming from the north. 

Assuming a wind speed of a mere 30 kilometres per hour, how can the wind 
direction affect the level of traffic noise I hear when the speed of sound is more 
than 1235 kilometres per hour? 

This apparent paradox, and its resolution, have been known since at least the time of 
Stokes [ij . The explanation given by Stokes is that this effect is produced by wind shear, the 
variability in the wind speed as a function of height. This gives rise to refraction, causing 



*DAMTP pre-print no. DAMTP-2009-75 
Email addresses: G.W.Gibbons@damtp.cam.ac.uk (G. W. Gibbons), C.M.Warnick@damtp.cam.ac.uk 
(C. M. Warnick) 
1 issue 2704 



Preprint submitted to Annals of Physics 



October 8, 2010 



sound rays to bend away from the ground in the upwind direction and towards the ground 
in the downwind direction. 

As a qualitative explanation of the process, this is perfectly satisfactory, see for example 
the discussion of Reynolds A quantitative discussion is given by Rayleigh [i[ however 
he assumes that the rays are normal to the wavefronts, which is not true in the presence of 
a wind. Rayleigh thus arrives at an incorrect equation for sound rays in a wind, as pointed 
out by later authors [3, |i[ . 

Modern approaches to this problem usually involve the use of numerical ray-tracing to 
plot the paths of sound rays [6]. The purpose of this paper is to provide a quantitative, 
analytic, discussion of this effect. We consider a stratified atmosphere whose sound speed 
and wind speed are allowed to vary with height. Making use of geometrical ideas discussed 
in [?J we show that provided the wind speed is small compared to the speed of sound, and 
that the sound speed does not vary rapidly with height, the rays are well approximated 
by trajectories of a charged particle moving in a uniform magnetic field on the hyperbolic 
plane. 

We use the equivalence between sound rays and charged particles to investigate two 
problems. One is that of traffic noise, outlined above. The other is the problem of the 
Mach envelope of a moving body in a stratified atmosphere with a small wind. We finally 
discuss how to lift the problem to that of finding null geodesies for a squashed anti-de Sitter 
spacetime, which allows us to exhibit the 50(3, 1) xl symmetries in a simple fashion. 



2. Fermat's principle for moving media 

Consider a disturbance u(x,t) obeying the wave equation 



d 2 Lij d d 



dt 2 dx i dxi 



u(x,t) = 0, (2.1) 



where h^(x) is a positive definite symmetric matrix. A sound wave propagating locally in 
the direction n l will have phase velocity h^n 1 ^ at x. Here hy is the matrix inverse of h % K 
For the case of an isotropic fluid, we have = c 2 (x)5 i: ' , with c the local speed of sound, 
but we will allow for the possibility that the speed of sound depends on the direction of 
propagation. 

Suppose we take a disturbance whose wavelength is short compared to all other length- 
scales. It is well known that the energy of the wave travels along rays a = at(r), which are 
extremals of the funtional 



T[a] = / v , //, / n' ; n / 'r/7 (2.2) 



Where r is any parameter along the curve. We note that the integrand is simply dt, so that 
rays obey Fermat's principle of least time, i.e. a sound packet will travel from P to Q along 
a path which minimisesE) the time taken among nearby paths, subject to the condition that 
the packet moves always at the local speed of sound. 



2 strictly extremises - the ray may in fact maximise the time taken 
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Now let us consider a disturbance in a fluid with local sound speed tensor h^(x) which 
is moving with a local velocity W l (x). This will obey a modified wave equation: 



0_ 

Of 



_d_ 

dx i 



h lJ 



d d 

dx % dxi 



u(x, t) = 0. 



(2.3) 



What is perhaps less well known is that the sound rays in this situation also obey a form of 
Fermat principle: a sound packet will travel from P to Q along a path which minimises the 
time taken among nearby paths, subject to the condition that it moves at the local speed of 
sound, relative to the flow. 

This means that a sound ray in a moving medium must solve a Zermelo navigation 
problem. This problem, proposed by Zermelo is to find the path between two points which 
minimises the time travelled, given that one moves at unit speed with respect to a 'wind' 
vector W. Taking the view that hij defines a metric, with respect to which sound moves at 
unit speed, the principle of least time tells us that a ray indeed solves a Zermelo problem. 

The general Zermelo problem of navigation on a manifold with metric h and wind W 
is considered in jsj. It is shown that this problem is isomorphic with finding the Finsler 
geodesies of a metric in the Randers class of Finsler metrics. Randers metrics take the form 



F : TM \ -> 



(x, y) i y ^a ij {x)y i yi + b { y l 



(2.4) 



Where a^jb 1 ^ < 1 is required in order that this be a good Finsler metric. The Randers 



metric whose geodesies solve the Zermelo problem with data {hij, W l 



Xhjj + WjWj 
X 2 



X ' 1 



hijW 3 , 



is given by 



hijW'WK 



(2.5) 



The condition a^W < 1 becomes hijW l W J < 1. This transformation is invertible and 
every Randers metric may be thought of as arising from a Zermelo problem. Geodesies of 
the Randers metric are (up to parameterization) the paths followed by a particle of unit 
charge, moving at unit speed in the magnetic field defined by F = db. For this reason, we 
will identify Randers one-forms which differ by an exact form b ~ b + dip for most of this 
paper, the exception coming in Section |5] where the parameterization will be relevant. 

The fact that the rays of (12. 3p obey Fermat's principle can be seen in two ways. In [9| 
it is shown that one may start with the Hamiltonian system following from the dispersion 
relation of (12.31) and show that the integral curves of the Hamiltonian obey the Euler- 
Lagrange equations of the time functional 



T\a] 



aija' l a'i + bia l ds 



(2.6) 



with a and b defined in terms of h and W by (12. 5p . Since the Hamiltonian is homogeneous 
of degree 1 in momenta this transformation is not the standard Legendre transformation. 
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Alternatively, in Fermat's principle is derived by considering (I2.3P as the d'Alembertian 
operator of a 4-dimensional Lorentzian manifold. Using the fact that the rays are determined 
by the principal part of the wave equation (which is conformally invariant) one may derive 
Fermat's principle by considering the null geodesies of this Lorentzian spacetime. In [?J the 
relation between the Zermelo, Randers and spacetime viewpoints was explored in greater 
depth. 

The link between magnetism and wave propagation in a moving background has also been 



explored in a different context by Berry et al. [10|. Here the relation is used to construct a 



water wave analogue for a quantum mechanical system. 
3. Traffic Noise 

We consider the following problem as a model for traffic noise in the vicinity of a mo- 
torway, with a stratified atmosphere such that the speed of sound varies with height above 
the ground and there is in addition a cross wind, also varying with height. We work in the 
upper half plane, with coordinates x, z and take the ground to be z = 0. We suppose the 
sound waves travel at a speed c(z) and there is a horizontal wind W = w(z)d/dx which 
vanishes at z = 0. Thus, we seek to solve the Zermelo problem with data 

The corresponding Randers data are 

_ dx 2 dz 2 _ w(z) 

fl %w(i-^)* + **) (i-=$y b= ~w^m dx - (3 - 2) 

If we make the assumptions that w(z)/c(z) <C 1 and c"(z)c(z)/c'(z) 2 <C 1, then the 
Gaussian curvature of the metric a is 

K —(c'(z) 2 + 2w'(z) 2 ). (3.3) 

We will therefore approximate the metric a in the neighbourhood of z = by a metric of 
constant negative curvature K = —[o 2 c + 2cr^), where we define for convenience a c = c'(0), 
a w = w'(0). The line z = may be shown to have geodesic curvature 

k g = a c (3.4) 

so in our approximation, the surface of the ground will be a curve of constant geodesic 
curvature in the hyperbolic plane. a c is positive when the speed of sound increases with 
height from ground level and negative if the speed of sound decreases with height. In the 
model proposed, the sign of a c determines on which side of the curve of constant geodesic 
curvature the atmosphere lies. 
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To first order in w/c, we find 



F = db = w'(z)n a (3.5) 

where \i a is the area form of the metric a. Thus in our approximation we will consider 
the sound rays to be the paths of a particle or unit charge, moving at unit speed on the 
hyperbolic plane with respect to a constant magnetic field a w . 

We use these observations to map the problem to one in the hyperbolic dislfl. This is 
the region \(\ < 1 of C, endowed with the Poincare metric 

ds 2 = V d(d(. (3.6) 

(i-ICI 2 ) 

This is a model of the hyperbolic plane, with constant curvature K = —1/p 2 . A particle of 
unit charge moving in a constant magnetic field /3, passing through the origin, moves along 
an arc of a circle of Euclidean radius 

R B = ± (3.7) 

Alternatively, we can characterise this arc as a curve of constant geodesic curvature (3. 

We will assume^ that the source of the noise is located at the centre of the disk. We then 
have the following building blocks for our geometric construction 

• Ground level corresponds to an arc of a circle of radius 

(3.8) 

through the origin. If a c > 0, so that the speed of sound increases with height, the 
atmosphere corresponds to the interior of this circle otherwise it is the exterior. By a 
rotation we may arrange that at z = the ground level is tangent to Imz = 0, with 
the atmosphere above the ground. 

Sound Rays are arcs of circles of radius 

(3.9) 

passing through the origin. If a w > 0, so that the wind blows from left to right above 
ground- level, the arcs should curve to the right, else they should curve to the left. 





3 In a later section, we will deal with the precise mappings which take our problem over to the hyperbolic 
disk, but we already have sufficient information to propose a geometric construction to determine received 
intensity as a function of distance from the source. 

4 this is easily arranged by making use of symmetries of the Poincare disk. 
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Figures [T] to HI show sample constructions for various values of a c and a w . The ground is 
shaded, and rays are shown emanating from the origin at equally spaced angles of 7r/12. Since 
the Poincare model is conformally flat, a source which radiates uniformly in all directions 
will radiate the same amount of energy between any two rays, 1/12 of the total energy 
radiated. We assume for simplicity that there is no significant transfer of energy between 
the bulk fluid motion and the sound waves. 




Figure 3: a c = 0, a w = .5 Figure 4: a c — — 1, a w = —2 



Figure [T] shows the case of vanishing wind, together with a sound speed which increases 
with height, such as may be caused by a 'temperature inversion'. In this case the rays are 
straight lines and all return to the ground at some point. Clearly one half of the power 
radiated returns to each side of the source. 
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Figure |2] shows the case where sound speed increases with height and there is also a wind 
shear, with the wind blowing from left to right. In this case, not all the radiated energy 
returns to earth. By counting the rays, we readily estimate that roughly 5/12 of the power 
is lost to the atmosphere, 5/12 is received downwind (i.e. to the right) and 2/12 upwind 
(left) of the source. 

Figure [3] shows the case where the sound speed does not vary significantly with height, 
but there is a wind blowing left to right. In this case, in the ray theory approximations, no 
sound is received to the left of the source, approximately 1/4 is received to the right and 
the rest is radiated upwards to the atmosphere. 

Figure S] shows the case where there is a decrease in sound speed with height, together 
with a strong wind shear, with the wind blowing from right to left. In this case there are 
competing effects between the changing sound speed, which would tend to refract waves 
upwards and the wind shear which would tend to bend rays downwards. We see that to the 
right of the source there is no sound received, while roughly 1/8 of the power is received to 
the left. In both Figure [3] and Figure @] there is a 'quiet zone', shown shaded in light grey, 
where no sound reaches the observer. 

We see then how a simple 'straightedge and compass' construction allows us to make 
quantitative predictions about the ratio of power transmission upwind and downwind in 
a shearing wind. In fact, with a little more geometry it is possible to calculate, in our 
approximation, the intensity of sound received as a function of distance from the source. 
Suppose a sound ray is emitted from the source at an angle 6 to the ground. It is a matter 
of simple circle geometry to calculate the point at which this ray again intersects the ground. 
Using the hyperbolic metric we can find the proper distance along the ground to this point, 
and we find that it is given by 




(3.10) 



where c = c(0). The received intensity is then simply proportional to (%) . We plot x 
against (||) 1 for the same values of <7 c ,a w as above in Figures [5] to [HJ 
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Figure 5: a c = 1, a w = 



Figure 6: cr c = 1, a., 



w 



= .5 



5 We would expect that for a full solution of the wave equation the sound field will not be zero but will 
in fact vanish exponentially in any 'silent' zone of the ray theory approximation. 
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Figure 9: a c = l,a w = 

In these figures, we have set cq = 1. The vertical scale is proportional to the intensity 
received. The plots agree with the rough ratios we found above for the energy upwind and 
downwind of the source. 

The approximations we have made are valid provided that a c x/co and <j w x/cq are small, 
so we can trust the plots in a neighbourhood of the origin. Qualitatively, we do not expect 
significant deviation from these plots, even outside the regime where the approximations are 
valid. 

We note that Randers metrics and Zermelo problems are in direct, one-to-one correspon- 
dence. We might therefore choose to take the view that the preceding calculations are exact, 
but for an atmosphere with a slightly different sound and wind speed profile. These profiles 
would match the assumed shape above in a region near the ground. 

4. Ray Plots 

In order to plot the paths of rays in the physical x, z coordinates, we need to explicitly 
construct the map to a region of the hyperbolic disk, of which we have made implicit use 
in the preceding section. We do this in the appendix below. In this section we use the 
mapping to pull back the ray paths from the hyperbolic disk to the original coordinates so 
as to exhibit the physical paths of the rays. Figures [9] to [12] show these paths for the same 
parameters as the previous section. 

We once again take c = 1. We note that the ray paths are broadly as one would expect 
from heuristic considerations [H, 0). In Figure [T2l we explicitly see that the rays terminate 
at z ~ .42. The horizontal line here is the pull back of the conformal boundary of H 2 , and 
the rays will take an infinitely long time to reach this height. Figures [10] and [IT] would 
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also exhibit such a phenomenon if the z-axis were extended. In practice of course, the 
approximations we have made are only valid in a strip about z = and would break down 
before the rays reached this height. 

One may verify roughly that these figures are consistent with the intensity plots given 
above. The energy absorbed by the ground between any adjacent rays is tt/12, which should 
be the area under the intensity plot between the same values of x. 



5. The sound field of a moving source 

We now change perspective somewhat and leave behind the problem of propagation of 
traffic noise to instead consider a moving body which radiates sound, for example an aircraft. 
In the classical case of a spatially homogeneous 2-dimensional atmosphere with no wind, a 
subsonic aircraft moving along the z-axis at constant speed has a soundfield that fills the 
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Figure 13: Example sound fields for a subsonic (t) and supersonic (b) body in a homogeneous atmosphere 
with no wind 



whole of IR 2 . For an aircraft moving at a supersonic speed, the disturbances are confined 
inside the 'Mach cone' and there is a zone of silence in front of the body. There is a 'Sonic 
Boom' at the boundary of this cone. These features are shown in Figure [T51 It is known that 
the introduction of a spatially varying sound speed gives rise to interesting, qualitatively 
different phenomena in the sound fields of fast moving objects 12]. We shall examine the 
further effects to be observed in the presence of a wind. 

In order to visualise the sound field, it is convenient to imagine that the object emits 
pulses of sound periodically, and to plot the loci of all these disturbances at some fixed 
time. We work in the approximation outlined in Section 121 where the wavelength of sound is 
assumed to be much shorter than other lengthscales in the problem and so the disturbances 
propagate along geodesies of the Randers metric F of (12. 4[ 12. 5p . The locus of a disturbance 
emitted a time t in the past is therefore a geodesic circle of radius t with respect to the 
Randers metric F, centred on the point at which the disturbance originated. 

We shall consider a body moving in the stratified atmosphere of Sections[3JH]at a constant 
speed parallel to the x-axis. After a Galileian boost, we may assume that the wind speed is 
zero at the height of the body, so that it moves along z = without loss of generality. In 
order to visualize the sound field as described above, we need to find the geodesic circles of 
the Randers metric defined by 

b = - a , ™ {Z \, , dx, (5.1) 
c 2 (z) — w 2 (z) 



c 2 (z) 



dx 2 dz 2 



centred on points of z = 0. Since d/dx is a Killing vector of the Randers metric, we may 
without loss of generality consider only those geodesic circles centred on the origin. 
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We once again assume that w(z)/c(z) <C 1 and c"(z)c(z)/c'(z) 2 < 1, so that the geodesies 
of the Randers metric are well approximated by the rays found in previous sections. We 
note that the time taken to traverse a curve 7(A) has two contributions: 

t = J dA(v^V¥ + 6<f) (5.2) 

but that (15. ip shows the bff contribution to be smaller than the first contribution by a 
factor of w/c. We may thus approximate the length of a curve with respect to the Randers 
structure a + b by the Riemannian length of the metric g£| We shall find the geodesic circles 
of radius t by moving a distance t, with respect to a, along the geodesies calculated in 
previous sections. We make use of the approximations above to map the problem to the 
hyperbolic disk, where the manifest rotational symmetries of the problem show that the 
geodesic circles we seek are in fact Euclidean circles centred on the origin. The relation 
between the Euclidean radius and the Randers radius is found by integrating a along an arc 
of a circle. Finally the transformations given explicitly in the Appendix are used to map the 
geodesic circle back to physical space. Combining this construction with the ^-translation 
invariance of the physical space, we may readily plot the sound field for various values of 
the parameters a c , a w and for differing speeds of the radiating source, expressed as a Mach 
number as v — Mc , with c the speed of sound at z = 0. 




Figure 14: Sound field for a body with M = 2 in an atmosphere with er c = 1, o~ w = 

Figures HH to [16] show some sample sound fields. We see that Figure HH is qualitatively 
very similar to that of [12(] which considers a moving source in a stratified atmosphere whose 
speed of sound varies as c(z) = (1 — z)~ x l 2 . This particular sound speed profile is chosen do 
that the rays are parabolae with vertical axes. For such an atmosphere, there is no 'zone 
of silence', essentially because rays which travel into the upper half plane and are refracted 
back downwards can 'overtake' the moving source. This is because the source is moving 
at a speed which is supersonic only with respect to the local sound speed, whereas in this 
atmosphere the speed of sound increases with height so that a ray travelling at a sufficient 



6 It might seem that we have removed any possible effect the wind might have by such an assumption. It 
should be noted however that the geodesies themselves depend on the wind. In this approximation, we keep 
the refractive effect of wind shear, while ignoring the transport effect of the wind, which will be negligible 
for reasonable wind speeds. 
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Figure 15: Sound field for a body with M = 2 (t) and M = 1.3 (b) in an atmosphere with <j c = 1, cr w = 0.5. 
Note the altered range of the lower figure. 




Figure 16: Sound field for a body with M = 2 in an atmosphere with a c — 0, a w = 0.5 

height will travel faster than the source. There is however a 'Mach envelope', so an observer 
may experience a sonic boom. The envelope has a cusp at the height at which the local 
sound speed equals the speed of the source. 

Figure [15] shows sound fields for a supersonic object moving in a stratified atmosphere 
with both a varying speed of sound and a wind. We find that there are two regimes. For 
a source moving much faster than the speed of sound, the field is similar to the Mach 
cone of a supersonic body in a homogeneous wind-free atmosphere, albeit distorted by the 
inhomogeneity. For a source which moves more slowly, but still above the local speed of 
sound, the picture is more similar to that of Figure (TU The reason for these different 
regimes is that the wind shear and varying speed of sound produce competing effects. The 
varying speed of sound tends to bend rays above the axis downwards, while the wind shear 
tends to bend rays travelling from right to left upwards. For a source travelling just over 
the local speed of sound, rays which travel into the upper half-plane and are refracted down 
may overtake the object as in Figure [HJ For a source travelling considerably faster than the 
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local speed of sound this cannot happen as any ray which moves far enough into the upper 
half-plane so that it moves faster than the source is refracted by the wind shear and never 
returns to z = 0. 

Figure [161 shows the case of an atmosphere whose sound speed is constant, but which has 
a varying wind. In this case, the sound field is similar to the Mach cone of the homogeneous 
wind-free case, but that the cone bends inwards slightly away from the axis due to the 
refraction caused by the wind shear. 



6. SL(2, M), squashed AdS3 spacetimes and integrability 

In [3] we explored a third vertex to the correspondence between the Zermelo problem and 
the problem of finding geodesies of a Randers metric. Both problems, and the correspondence 
between them, fit naturally into the problem of finding the null geodesies of a conformally 
stationary Lorentzian spacetime. Firstly we note that any equivalence class of metrics [g] 
which are conformally stationary locally have a representative of the form 

gn = — (dt — bi(x)dx 1 ) 2 + dij(x)dx' 'dx 1 (6.1) 

with [gn] = [g], which we call the Randers form of the metric. If we choose to parameterize 
the null geodesies of [g] by t, then they are in fact geodesies of the Randers metric 

F(x, y) = y/oijixWyi + hy 1 (6.2) 

parameterized by unit Finslerian length. This choice of representative naturally picks out a 
Riemannian metric a together with a one-form b. 

Another way to write a manifestly stationary spacetime, often referred to as the Painleve- 
Gullstrand form, is given by 

gz = -dt 2 + hij(x) (dx i + W\x)dt) (dx j + W j (x)dt) . (6.3) 

We show in [?J that [gz] = [gn] if and only if the Randers metric (I6.2p solves the Zermelo 
problem with data (hij, W % ). The condition that the Randers metric be strongly convex is 
also necessary to ensure that the conformal factor is nowhere singular. We thus see that the 
Randers and Zermelo pictures arise as different ways to write a conformally stationary class 
of conformal metrics. 

We now have a prescription to lift the Finslerian geometry of the Randers structure to 
the more familiar Lorentzian geometry. The problem we have explored in depth above is 
that of a constant magnetic field ft on the hyperbolic plane with curvature — k 2 . Following 
Q, we note that the lifted metric in the Randers form may be written 

9R = o ( 6 - 4 ) 
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where the pi are right invariant one-forms on SL(2,. 

i ndt 1 + r 2 . 
P = -1T + 



2 



P 1 -r 
2 

p 2 = (r cos(nt / f3) d(p — sm(nt / f3) dr) 

1 — r 2 

2 

p 3 = (r sin( nt / P) dfi + cos( nt / P) dr) . (6.5) 

1 — r 2 



These arise by parameterizing U G SL(2, M) as 



U 



^ W+X V -Y X 



■V-Y W-X 



X + iY = J r _ e »(^-W)/(2/8) ; 



(6.6) 



i v I vv — y\ j W + %V — ^ 1 _ r 



— < 



The standard procedure of writing the Maurer-Cartan form 

df/f/- 1 = p V i; (6.7) 

with 7j a basis for s/(2, R) given in terms of the standard Pauli matrices as {f 02, |o"3, — |cri} 
yields f E3|) . 

Thus the sound rays we found are projections of null geodesies of a right-invariant metric 
on SL(2,M). The ease with which we were able to construct the rays is then seen to be 
a consequence of symmetry. The metric (16. 4p certainly admits all the left-invariant vector 
fields of SL(2,M) as Killing fields. These are readily constructed as the dual basis to the 
left-invariant one-forms which satisfy U~ 1 dU = AVj and are given by: 

L - 9 

1 + r 2 d 1 — r 2 / d d 

L 2 = — — cos 0— H — r sin <p— - p cos 0— 

2r o<p 2r \ or at 

1 + r 2 d 1-r 2 ( d d\ 

In addition, since the coefficients of p 2 ,p 3 in the metric are the same, there is a further 
Killing field 

R x = P^, (6.9) 

so that the full symmetry group of the metric is (SL(2, M) x M)/Z2. The Z2 quotient cannot 
be seen at the level of the Lie algebra, but one sees upon exponentiating the algebra that it 
is in fact a double covering of the symmetry group. Using the isomorphism 

50(2,1) = SL(2 1 R)/Z 2 (6.10) 
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we can also identify the symmetry group as 5*0(2, 1) x R with the S0(2, 1) group rep- 
resenting the symmetries of the hyperbolic plane and the factor IR the time translation 
symmetry. There are many similarities between the situation we are dealing with and that 
of right-invariant metrics on S 3 = SU(2), which should come as no surprise since the groups 
are closely related by analytic continuation. The metrics we have exhibited are a natural 
analogue of the biaxial Berger metrics on S 3 , with symmetry group (SU(2) x [/(1))/Z 2 . 
The metrics (16 .4p are often referred to as describing a squashed anti-de Sitter geometry 



1 31 ] . For (3 2 > k 2 we have a family of metrics similar to (a factor of) the Godel universe, 
which include the Godel universal as the special case (3 2 = 2k 2 and which all admit closed 
timelike curves (CTCs). For (3 2 < k 2 the spacetimes do not admit CTCs. 

After approximation, we found above that the problem of finding sound rays for a strat- 
ified atmosphere with a wind shear is equivalent to that of finding the motion of a charge 
in a constant magnetic field of (3 = a w in a hyperbolic space with k 2 = a 2 + 2a^ u) with 
<j c = c'(0), <j w = w'(0). Clearly we will always be in the regime where /3 < n so do not 
expect any CTCs in the lifted spacetime. This is related to the fact that any two points in 
the space may be joined by a Randers geodesic and these properties are explored in more 
detail in 



7. Conclusion 

We have examined the old problem of ray tracing for sound waves in a wind. We find that 
in a certain physically plausible approximation the ray paths are mapped to the trajectories 
of a charged particle moving in a uniform magnetic field on the hyperbolic plane. We have 
exploited this fact to produce plots of the intensity of sound received from a point source at 
ground level in the presence of a wind and varying speed of sound. We have also constructed 
the approximate Mach envelopes for a moving body in the presence of a wind. We have 
discussed the symmetries of this model and how they relate to the problem of finding null 
geodesies of a family of right invariant metrics on SL(2,M.) which are often referred to as 
squashed AdS metrics. 
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Appendix - Mapping to H 2 . 

We start with the metric a from above: 

dx 2 dz 2 



7 strictly speaking the Godel universe is the product of this space with a fiat direction 
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and make a near-identity chang^l of the vertical coordinate 



dy= x \l-^dz (.2) 



so that the metric may be written in the form 

dx 2 + dy 2 dx 2 + dy 2 



3) 



c%)(l-^?) P(V) ' 

At this point we will explicitly approximate f(y) as 

f(y) = acos(/3y - 7) (.4) 

which for 

is accurate to 0(y 3 ). After a shift in the y coordinate defined by 



y = y-jj (-6) 



we have 

dx 2 + dy 2 d(d( 



- ( 7) 

a 2 cos 2 /3y a 2 cos 2 [/?(( - C)/2i] 

where we define £ = x + iy. This is a (perhaps unusual) metric for the hyperbolic plane of 
curvature —a 2 (3 2 , where the strip — 7r/(2/3) < y < n/(2f3) is a complete copy of the plane. 
We map this strip onto the standard unit disk with the conformal mapping 

tu = tanh — (.8) 

so that the metric a becomes 

4 dudu 

a = ^(i- M 2 ) 2 [) 

The ground is mapped to a curve of constant geodesic curvature passing through oj — ±1. 
A final Mobius transformation leaves the metric unchanged, but shifts the ground so that it 
passes through the origin 



, pco + i 



where p = V2— + , I + 2^f. (.10) 



— no + p a c y a 



8 in practice we may ignore this change of variables as y = z to 0(y 3 ), which is the order of the approxi- 
mations we make 
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